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We study the effect of a torque induced by the local molecular field on the phase ordering dynamics 
of the Heisenberg model when the total magnetization is conserved. The torque drives the zero- 
temperature ordering dynamics to a new fixed point, characterized by exponents z = 2 and A ~ 5. 
This 'torque-driven' fixed point is approached at times such that g^t ^ 1, where g is the strength 
of the torque. All physical quantities, like the domain size L{t) and the equal and unequal time 
correlation functions obey a crossover scaling form over the entire range of g. An attempt to 
understand this crossover behavior from the approximate Gaussian Closure Scheme fails completely, 
implying that the dynamics at late times cannot be understood from the dynamics of defects alone. 
We provide convincing arguments that the spin configurations can be decomposed in terms of defects 
and spin- waves which interact with each other even at late times. In the absence of the torque term, 
the spin waves decay faster, but even so we find that the Gaussian closure scheme is inconsistent. 
In the latter case the inconsistency may be remedied by including corrections to a simple gaussian 
distribution. For completeness we include a discussion of the ordering dynamics at Tc, where the 
torque is shown to be relevant, with exponents z = A — e/2 and X = d (where e = 6 ~ d). We show 
to all orders in perturbation theory that A = d as a consequence of the conservation law. 
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I. INTRODUCTION 

When a many-body system like a magnet or a binary 
fluid is quenched from its disordered high temperature 
phase to its ordered configuration at low temperatures, 
the slow annealing of "defects" (interfaces in binary flu- 
ids, vortices (hedgehogs) in XY (Heisenberg) magnets) 
separating competing domains, makes the dynamics very 
slow. The system organizes itself into a self similar spatial 
distribution of domains characterized by a single diverg- 
ing length scale which typically grows algebraically in 
time L{t) ^ This spatial distribution of domains is 
reflected in the scaling behavior of the equal-time corre- 
lation function C(r, t) ~ f{r/L{t)). The autocorrelation 
function, A{t) ~ L{t)^^ is a measure of the memory of 
the initial configurations. The exponents z and A and the 
scaling function f{x) characterize the dynamical univer- 
sality classes at the zero temperature fixed point (ZFP) 

0. 

The above phenomenology suggests that the asymp- 
totic dynamics of the order parameter is dominated by 
the dynamics of its defects, and that bulk fluctuations 
(concentration waves in a binary fluid, spin waves in a 
magnet) relax fast and decouple from the dynamics of 
defects at late times. This picture is at the heart of re- 
cent approximate theories such as the Gaussian Closure 
Scheme 

But is this picture accurate ? In this paper we shall 
study the very realistic example of the conservative dy- 
namics of a Heisenberg magnet driven by a torque in- 
duced by the local molecular field, and show that the 
longer-lived spin waves couple to the defects even at late 



times, driving the system to a new fixed point. The new 
'torque-driven' fixed point characterized by z = 2 and 
A ~ 5.05, is accessed after a crossover time tc ^ l/g^ 
(where g is the strength of the torque). Crossover scal- 
ing forms describe physical quantities at late times (like 
the domain size L{t, g) and correlation functions C(r, t, g) 
and A(t, g)) for all values of g. In the absence of the 
torque, the spin- waves decay faster, but even so we find 
that the Gaussian Closure Scheme is internally incon- 
sistent. This inconsistency may however be rectified by 
including leading corrections to the gaussian distribution 
(as suggested by Mazenko for the dynamics of the 
conserved scalar (Ising) order parameter). 

For completeness we also study the effects of including 
the torque in the dynamics following a quench to the 
critical point T^- As reported in earlier studies [Q, the 
torque is relevant with exponents z — A — e /2 and X = d 
(where s = 6 — d). We show to all orders in perturbation 
theory that X = d which follows as a consequence of the 
conservation of total magnetization |^-^ . 



II. HEISENBERG MAGNET AND 
PRECESSIONAL DYNAMICS 

The order parameter (f> (whose components are 0^ with 
a = 1, 2, 3) describing a coarse-grained spin density in a 
Heisenberg ferromagnet in three dimensions experiences 
a torque from the joint action of the external field (if 
present) and the local molecular field. In response the 
spins precess with a Larmour frequency f^L about the 
total magnetic field. Coupling to various faster degrees 
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of freedom like lattice vibrations or electrons, causes a 
dissipation in energy and an eventual relaxation towards 
equilibrium. 

This dynamics follows from the generalized Langevin 
equation and the Poisson algebra S , 



dt S(f>a 



SF 



(1) 



The noise ff arising from the heat bath has mean zero 
and is conservative, 

(77„(x, i)7,^(x', t')) = 2 ksTTd^p V25(x - x')(5(t - t') . 

(2) 

The free-energy functional F is taken to be of the 
Landau-Ginzburg form, 



Fl 



= d 



• (3) 



The second term in Eq. is clearly the torque Af x H, 
where H = ~SF/6<j> is the local molecular field. 

Both the inertial term (by virtue of F being rotation- 
ally invariant in spin space) and the dissipation conserve 
the total spin, and so the full equations of motion (|l|) 
also conserve the total spin. 

Since the noise correlator is proportional to tempera- 
ture, we may drop it in our discussion of zero temperature 
quenches. We then scale space x, time t and the order 
parameter as 



X, t 



TrH 



G a 

to obtain the equation of motion in dimensionless form, 



dt 



(0X V2^) 



(4) 



The dimensionless parameter g = (r2icr/r)(ru) is 
the ratio of the precession frequency to the relaxation 
rate. Setting ~ 10^ Hz and T ~ 10*^ - 10^° Hz, gives 
g in the range of ^ 10^^ — 10. 



III. PHASE ORDERING DYNAMICS AT T = 

Let us now prepare the system initially in the para- 
magnetic phase and quench to zero temperature. We 
study the time evolution of the spin configurations as 
they evolve according to Eq. (||). We calculate the equal 
time correlator, 



C(r,i) = (<^(x,i)-(^(x + r,t)) 
and the autocorrelator. 



(5) 



C(0, h = QM=t) = A{t) = (0(r, 0) • 0(r, t)) , (6) 

where the angular brackets are averages over the random 
initial conditions and space. At late times these correla- 
tors should attain their scaling forms 



C(r,t)~/(r/L(i)) 



and 



A{t) ^ L{t) 



-A 



(7) 



(8) 



The length scale L{t), which is a measure of the distance 
between defects, may be evaluated either from the first 
zero of C(r, t) or from the scaling of the energy density, 
£ = y J dr {{ V0(r, i) )^ ) ^ L{t)^'^ , and grows with time 
as L{t) ~ t^/^ . We compute the scaling function f{x), the 
dynamical exponent z and the autocorrelation exponent 
A by simulating the Langevin Eq. (0). 



A. Langevin Simulation 

The Langevin simulation is performed by discretizing 
Eq. on a simple cubic lattice (with size TV ranging 
from SO'^ to 60^) and adopting an Euler scheme for the 
derivatives The space and time intervals have been 
chosen to be Aa; = 2.5 and Ai = 0.2. With this choice 
of parameters, we have checked that the resulting cou- 
pled map does not lead to any instability. We have also 
checked that the results remain unchanged on slight vari- 
ations of Ax and At. Throughout our simulation we have 
used periodic boundary conditions. 

The correlation functions Eqs. (||), (H) are calculated 
for values of g ranging from to 1. Measured quantities 
are averaged over 5—10 initial configurations. The initial 
configurations are taken from two ensembles, both in the 
disordered phase. In ensemble A, (t){t — 0) is uniformly 
distributed within the volume of a unit sphere centered 
at the origin, (j) at different spatial points are uncorre- 
lated. In ensemble B, 0(t = 0) is uniformly distributed 
on the surface of a unit sphere centered at the origin. 
(j) at different spatial points are again uncorrelated. We 
consider these two initial conditions to check if the late 
time dynamics is insensitive to the choice of initial en- 
semble (as long as they do not introduce any long-range 
correlations). 

We first report simulation results for ensemble A. 

Figure 1 is a scaling plot of C(r, t) versus r/L{t) for var- 
ious values of the parameter where Lit) is extracted 
from the first zero of C(r, t). Note that the scaling func- 
tion for g = Q is very different from those for g > 
further the g > scaling functions do not seem to de- 
pend on the value of g. This suggests that the dynamics 
crosses over to a new 'torque-driven' ZFP. This is also 
revealed in the values of the dynamical exponent z. In 
Fig. 2, a plot of Lit) versus t gives the expected value of 
z = A when g = Q. For g > 0, we see a distinct crossover 
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from z = 4 when t < tc{g) to z = 2 when t > tc{g)- The 
crossover time ic(ff) decreases with increasing g. The 
same z exponent and crossover are obtained from the 
scahng behaviour of the energy density e. 

To make sure that our resuhs are not affected by finite 
size, we compute 3 relevant time scales (shown in Table 
1 below) — (1) tc{g), the crossover time from a t^^'^ to 
a t^/^ growth, (2) ts{g), the time at which asymptotic 
scaling begins, (3) tjs, the time at which finite size effects 
become prominent. It is clear from the Table that tc < 
tg < tfs, as it should be if our data is to be free of finite 
size artifacts. A general rule-of-thumb is that finite size 
effects start becoming prominent when the domain size 
gets to be of order 1/3 the system size, and we see from 
Table 1 that L^ax/N is comfortably less than 1/3. 



Tabic 1 
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FIG. 2. Log-Log plot of L{t). At g = (o) we find 
that z ~ 4 (line of slope 0.25 drawn at the bottom for 
comparison). At g ^ {g = 0.1(*), 0.3(0), 0.5(A)), z 
crosses over from 4 to 2 (line of slope 0.5 drawn at the 
top). 
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FIG. 1. Scahng plot of C(r,t) for N = 50^. The scal- 
ing function f{x) changes as g is varied from (7 = (o) to 
.g^0(g = 0.1(*), 0.3(A), 0.5(n)). 



The last column in Table 1 shows fmin, the value of the 
scaling function evaluated at the first minimum as a func- 
tion of g. It is easy to see why fmin{g) < fminig = 0), 
since the precession of the spins about the local molecu- 
lar field would cause spins from neighboring "domains" 
to be less anti-correlated. This is borne out by comput- 
ing the spin-wave correction to an approximate form of 
C{r,t;g = 0) (given in Eq. (p9|), more on this later) to 
quadratic order in the spin- wave amplitude 

The autocorrelation function A{t) is calculated for 
g = 0,0.2 and 0.3 (Fig. 3). The simulations have been 
done on a lattice of size 60'^ and averaged over 10 initial 
configurations (we have to average over a large number 
of initial configurations for smoother data). The A expo- 
nent extracted from the asymptotic decay of A{t) clearly 
suggests a crossover from A = 2.2 to A « 5.05. The nu- 
merical determination of A is subject to large errors [p^P] 
and is very sensitive to finite size effects, and so we have 
to go to very late times and hence large system sizes to 
obtain accurate results. 
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FIG. 3. Log-Log plot of A{t) vs t for g = 
0(o), 0.2(+), 0.3(n). Solid line on top has the form 
a/t^^^ where A = 2.19 and z = 4 (corresponding to the 
5 = fixed point) while the one below has a A = 5.05 and 
z = 2 (corresponding to the 'torque-driven' fixed point). 



To make sure that we collect asymptotic data un- 
tainted by finite size, we compute two time scales (Ta- 
ble 2) — (i) tfit{g), the time beyond which A{t) can be 
fit with a power law a{t + fo) (ii) i/s, the time at 
which finite size effects on A{t) become prominent. The 
crossover time tc was displayed in Table 1. 

To determine t we plot an effective exponent Ag/ / — 
—td(\og A{t)) / dt as a function of 1/t. The derivative is 
calculated numerically with a (5t = 15 (in units of the 
time discretisation At). We see from Fig. 4, that at late 
times t > tfs, Xeff crosses over to being a decreasing 
function of time, clearly a finite size effect. This estimate 
of tfs is not very sensitive to the choice of St changing 
by 1% (for g = 0.2) and 3.5% (for g = 0.3) as 6t changes 
by 5 units. Note that finite size effects in A{t) appear 
earlier than in C{r,t). 

It is seen from Table 2 that tfu < tfs, as it should if 
we are to have an accurate determination of A. 

Table 2 
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FIG. 4. Xeff versus 1/t for g = O.O(o), 0.2(*), 0.3(+). 
Finite size effects set in when Ae/ / starts becoming a de- 
creasing function of time. For 5 = we do not see any 
finite size effects in A within our simulation times. 

The last column of Table 2 lists the value of A as a func- 
tion of g. The data presented and the plot in Fig. 3 
clearly support a crossover from A = 2.2 at 5 = to 
A = 5.05 at g ^ 0. The values of A satisfy the bound 
derived in | |lO|] . 

We now present results of the Langevin simulation for 
initial conditions taken from ensemble B. We find that 
the value of z, the form of the scaling functions f{x) (Fig. 
5) and the decay of the autocorrelation function A{t) 
(Fig. 6) are insensitive to the choice of initial conditions. 
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FIG. 5. Scaling function f{x) vs x for g = 0.3 using 
ensembles A(o) and B(A). 
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Since the initial condition B sets the magnitude of the 
spins to its T = equihbrium value, the crossover time 
tc is smaller than for ensemble A. For the same reason 
the domain sizes computed using ensemble B are larger 
than that of A. 
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FIG. 6. Log-Log plot of A{t) for g = 0.3 using ensem- 
bles A(o) and B(-l-). A power law a/i^/^ with A = 5.05 
is displayed for comparison. 



where Mq, a and are the equilibrium magnetization, 
surface tension and spin mobility respectively. Beyond a 

crossover time given by (5) ^ (L/Afoi^)^ ^ ^/g^i simple 
dimension counting shows that the dynamics crosses over 
from z = 4 to z = 2 in conformity with our numerical 
simulations. 
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FIG. 7. Scaling plot of ?/ = L{t,g)/t'^/^ versus x = tg'^ 
for g = 0.03(o), 0.05(0), 0.07(A), 0.09(o), 0.10(*). The 
solid line of slope 0.25 is the theoretical estimate of the 
asymptotic form of the scaling function as a; ^ 00 (see 
text). 



B. Crossover Phenomenon 

It is clear from the last section, that though the asymp- 
totic dynamics is governed by the new 'torque-driven' 
fixed point, the dynamics at earlier times t < follows 
the 5 = behavior. This suggests that the dynamics 
for arbitrary may be analyzed as a crossover from the 
(7 = fixed point characterized by (z = 4, A ~ 2) to the 
torque-driven fixed point where (z = 2, A w 5). 

A simple scaling argument encourages us to think of 
such a crossover scenario. On restoring appropriate di- 
mensions, the dynamical equation Eq.(0) can be rewrit- 
ten as a continuity equation, 

d${v,t)/dt^~v ■] (9) 

where the "spin current" is 

Ja - -r (^V^ + ^e„^^,^^V0^ j . (10) 

From a dimensional analysis where we replace ja by the 
'velocity' dL/dt, we find 



The crossover physics is best highlighted by numeri- 
cally demonstrating crossover scaling of the domain size 
L(t, g) and the correlation functions C(r, 17) and A{t^ g). 

For instance, Eq. ( pT| ) suggests that the domain size 
obeys the scaling form L{t,g) — t^/^Sm{tg'^) where the 
crossover function Sm{x) is determined from the tran- 
scendental equation, 

x^''^s^{x) - ln(l -f x^/^s^) - 2x = . (12) 

We shall now argue (and then confirm numerically) 
that the above scaling form holds in general. Scaling 
r r/b, t t/b^ and g ^ g/b^^, scales the domain size 

by 

L{t,g) = bs{t/b^g/by^) (13) 

where yg is the scaling dimension of g. We choose b such 
that t/b^ — 1, which implies 

L{t,g)^t'^^.sig/ty^^^). (14) 

Setting g = gives L{t, g = 0) — t^^^s{0), telling us that 
z = 4. Thus the scaling form Eq. (^^ is governed by the 
5 = fixed point. We therefore need to evaluate yg at 
this 5 = fixed point. We determine yg by noting the g 
contribution to Eq. (^ 
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X /i 
2 



(15) 



where the last relation is obtained by demanding lo- 
cal equilibrium (Gibbs-Thomson) on the chemical po- 
tential fl. Thus equating dimensions, [g\ = [i~^][L^] = 
= [L^^] leading to j/g — —2. The crossover scal- 
ing form for the domain size can now be read out from 
Eq. (0), 



L{t,g)^t^'\s{gh). 



(16) 



The X ^ oo asymptote of s(x) can be obtained by de- 
manding that we recover the 'torque-driven' fixed point 
behavior, which forces s{x — > oo) ~ x^^^. 

We will now check whether this crossover scaling form 
is seen in our Langevin simulation. If the above proposal 
is true, then the data should collapse onto the scaling 
curve s{x) when plotted as L{t, g)/t^^'^ versus tg^. Fig- 
ure 7 shows the results of the numerical simulation — 
the data collapse is not good away from the asymptotic 
regimes. To see a better data collapse away from either 
fixed point, it is necessary to include corrections to scal- 
ing. 

Corrections to scaling come from two sources — (i) 
finite time effects and (ii) nonlinear corrections to the 
scaling fields ||ll|. Finite time corrections can be incor- 
porated by introducing finite-time shift factors t t — to, 
which can be neglected in the t —i- oo limit. Nonlinear 
corrections to scaling are incorporated by constructing a 
nonlinear, analytic function g(g) of the physical fields g, 
such that it reduces to g in the limit g ^ 0. The simplest 
choice of such a function is 



K9) = 



eg- 



1 + cg^' 

leading to a nonlinear scaling variable 
x={gig)fit-to). 



(17) 



(18) 



The data plotted with respect to this nonlinear scaling 
variable shows a much better collapse (Fig. 8) when c 
is chosen to be around —1.5 (in the Figs. 8 - 10, the 
finite time shift was taken to be 0). The simple mean- 
field estimate Sm{x) plotted for comparison ( Eq. (p^), 
is exact only at the asymptotes. 



y 



-0.24 - 
-1.5 



-0.5 0.5 1 1.5 2 



FIG. 8. Plot oiy = L{t,g)/{t - to)^/'^ versus x when 
cfn —1.5. The point o on the y axis, represents the value 
of ?; as x — > 0. 

We have seen in the last section that the equal time cor- 
relation function C (r, t, g) is unaltered when scaled with 
the domain size L, and so we expect it to have the fol- 
lowing scaling behavior 



C{r,t,g)^f{r/L,t/L^g/Ly^), 



(19) 



where z is the dynamical exponent at the g = fixed 
point and yg is the scaling dimension of g. L is the size 
of the domain, given by Eq^ 10) . This readily leads to a 
two variable scaling form |l^ , 



C(r,i,g) = /(^,tg2) 



(20) 



with scaling variables p — rjt^l^ and x — tg^. When 
X = and x ^ oo then f{p,x) = fo{p) and f{p,x) = 
frip) respectively, where foip), frip) are the asymptotic 
scaling functions at 5 = and .g ^ 0. Again in terms of 
the nonlinear scaling variables x and p = r /(t — io)^^^, 
we find a very good collapse of the data for c « —1.2 
(Fig. 9). 
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FIG. 10. Log-Log plot oiy = A{t,g)/t-^°/'^ versus x 
for g — O.l(o), 0.2(-|-), 0.3(n) showing data collapse for 
c ~ —1.1. The scaling function aysmptotes to a line of 
slope Ao/4 — At/2 = —1.95 as i ^ cxa. 

The above discussion clearly indicates that for times 
t ^ tc{g) ~ 1/5^, the dynamics is affected by the 5 = 
fixed point while for t 3> tc{g) ~ 1/ff^, it follows the 
'torque- driven' fixed point. Our scaling analysis suggests 
the following renormalization group flow diagram, 



x = 
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C. Failure of Mazenko Closure Scheme : Interaction 
of Defects with Spin Waves 



FIG. 9. C{r,t,g) versus i at p = 0.50, 0.82, and 1.50 
for g = 0.03(o), 0.05(o), 0.07(0), 0.09(A), O.l(-I-), 0.3(*) 
showing data collapse for c « —1.2. 

Similar arguments suggest that the autocorrelation func- 
tion satisfies the scaling form 

A{t,g)^t-^»^^a{tg^), (21) 

where a{x = 0) = oq is a constant, and Aq ~ 2.2 is 
the value of the autocorrelation exponent at 5 = 0. As 
X — > 00, the scaling function a{x) should asymptote to 
a{x) ~ j,>^o/4:~^t/2 ^ where At ~ 5.05 is the exponent 
at the 'torque-driven' fixed point. This expectation is 
borne out by the numerical simulation (Fig. 10), where 
we have again used the nonlinear scaling variable x for 
better collapse. 
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We would like to know if the crossover phenomenon 
described in the last section can be understood from 
certain approximate theories of phase ordering of con- 
served vector order parameters. In particular could we 
use such theories to calculate the crossover scaling func- 
tions and the correlation functions at the 'torque-driven' 
fixed point. The Gaussian Closure Scheme introduced by 
Mazenko |^ has been considered a very successful theory 
to compute scaling functions of conserved vector order 
parameters, and it is to this we turn our attention. 

The method consists of trading the order parameter 
(/)(r, t) which is singular at defect sites, for an everywhere 
smooth field m(r, i), defined by a nonlinear transforma- 
tion, 

<^(r,t) = a(77i(r,t)) . (22) 

The choice for the nonlinear function a is dictated by the 
expectation that at late times, the magnitude of (j) satu- 
rates to its equilibrium value almost everywhere except 
near the defect cores. This suggests that the appropriate 
choice for a is an equilibrium defect profile, 

^\/la{rn{r,t))^V'{a{m{r,t))) , (23) 

where V'{x) = —x + {x ■ x)x. The auxiliary field rfi 
now has the natural interpretation as the position vector 
from the nearest defect core. Implicit in this choice is 
that smooth configurations such as spin waves relax fast 
and so decouple from defects at late times. The simplest 
nontrivial solution of Eq. ( |23| ) is the hedgehog configura- 
tion, 

'^(™(^'^» = |1^||^(H), (24) 
where 17(0) = and g{oo) = 1. 
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Equation (^) can be used to derive an equation for 
the correlation function C(12) = ((/)(ri,ti) • (/)(r2,t2))- 
Substituting for $ (Eqs. (|2|), (|2|)) in the right hand 
side of the resulting equation, we get 

dtC{12) = -Vl [VlC{l2) - { a{m{2)) ■ V'{a{rn{l))) )] 

+ g{a{rn{2)) ■ a{m{l)) x\/la{rn{l))) . (25) 

The Gaussian Closure Scheme assumes that each compo- 
nent of m(r, t) is an independent gaussian field with zero 
mean at all times. This implies that the joint probability 
distribution P(12) = P(m(l),m(2)) is a product of sep- 
arate distributions for each component and is given by 



We will first question the Gaussian assumption, on the 
lines suggested by Yeung et. al. Q in the case of a con- 
served scalar (Ising) order parameter. We will do this for 
the case when g = 0, the g ^ analysis follows similarly. 

The equal time correlation function may be derived 
from Eqs. @ and takes the form U 



Ztt 



1 2 



F[-- ^-r 

^ 2' 2' 2'^ 



(29) 



where B{x, y) and F(a, 6, c; z) are the Beta and hyperge- 
ometric functions respectively and 7 is given in Eq. ( p7| ) . 
We may expand the hypergeometric function as a power 
series in 7 urn and then take its fourier transform, 



exp 



where 



and 



, <(2) 
2(l-72)V^o(l) ^o(2) 

27TOa(l)mQ.(2) 

' VSo{l)So{2) 



27rV(l-72)5o(l)5o(2) 



(26) 



7 = 7(12) 



y/So{l)So{2) 



(27) 



The joint distribution has been written in terms of 
the second moments 5o(l) = and Co(12) = 

(m„(l)m„(2)). 

With this assumption, the right hand side of Eq. (|2£ 
simplifies to. 



dC{12) 
dti 



-V 



V^C{12) 



7 dC{l2) 



25o(l) dj 
g { ct(to(2)) • ct(to(1)) X V^a{m{l)) ] 



(28) 



where the laplacian is taken with respect to ri. With 
the joint probability distribution given by Eq. (|2^), it is 
clear that the last term in the above equation vanishes, 
implying that the torque is irrelevant at late times. This 
result of the Gaussian Closure Scheme, is in direct con- 
tradiction to the results of the last two sections. 

What has gone wrong ? There are two possible sources 
of error 

1. The Gaussian assumption for the probability dis- 
tribution of m is invalid. 

We show below that while the Gaussian assump- 
tion leads to an internal inconsistency, it may be 
remedied by considering corrections to the gaus- 
sian distribution. This however does not solve the 
above contradiction. 

2. The order parameter (j) cannot be written in terms 
of the defect field m alone. 



S(k, t) 



p=0 



dki . . . dk2p+i 



ap7ki(t)7k2(i)---7k2p+i(i) 



5(k + ki 



k2p+i) 



(30) 



where the spectral density 7k is the fourier transform of 
7(r, t) and the expansion coefficients. 



9 [r(p + i/2)]^ 

87r3/2 r(p + 5/2)p! 



B 2 



(31) 



are strictly positive for p > 0. If Eq.(pOD has to satisfy 
the conservation law S{k — 0,i) — 0, it is clear that 
7k (i) should be negative at some values of k. This is 
inconsistent with the definition Eq. (^^ which implies 
7k (i) > for all k. This definition is a consequence of 
the Gaussian approximation. 




FIG. 11. The spectral density 7(fc,t) at t = 3600 be- 
comes negative for < k/km < 0.5 and for 1.5 < k/km < 
3.0 (inset). 

To determine the range of values of k for which 7k 
is negative, we numerically evaluate the fourier trans- 
form of 7(r, t) after inverting Eq. (p^). This is prone to 
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numerical errors because of statistical errors in our com- 
puted C(r, i). For instance, a numerical integration of 
/ dr C{r, t) gives a nonzero value whereas it should be 
identically zero because of the conservation law. This is 
reflected in large errors in 7(k, t) at small k. We therefore 
adopt the following procedure. We fit a function Cf{x) 
to the equal time correlation function C{r, t) and use this 
to extract 7(k, t) from the Eq. (29). The fitting function 
has been taken to be 



Cf{x) 



sin(a;/L) 



1 



' a; \ 2 



exp[-fe(a;/L)2] (32) 



which is similar to the analytic form given in Ref. p5[ . 
Note that only b and L are independent fitting param- 
eters, a is determined from the condition Sf{k = 0) = 
0. This function with L = 1.5106 ± 1.01 x 10"'' and 
h = 0.0202 ± 2.14 X 10--* gives a very good fit to Cir, t) 
upto the fourth zero of the function. We observe (Fig. 11) 
that the spectral density, which should be a strictly pos- 
itive function of its arguments, becomes negative for 
k/km < 0.5 {'y{k,t) is peaked at km) and in the range 
1.5 < k/km < 3.0. 

Our demonstration suggests that a purely gaussian 
theory for the distribution of m is internally inconsis- 
tent. This may however be remedied by considering cor- 
rections to the purely gaussian distribution, as suggested 
by Mazenko for the scalar (Ising) order parameter. 

In order to help us understand the nature of the cor- 
rections, let us first numerically evaluate the probability 
distribution of rh. We determine m by choosing (7(|m|) in 
such a way as to make Eq. (|2j) invertible. A convenient 
choice is 



(j) = d{rh) 



(33) 



We now compute the asymptotic single point probability 
density P(mi(r, t)) on a 50*^ lattice averaged over 18 ini- 
tial configurations for g = 0, 0.3, 0.4 and 0.5. The prob- 
ability density obeys a scaling form at late times (Figs. 
12 and 13), P(mi,t) = P{mi/L{t)), where the length 
scale L{t) = \/ (m^) ~ i)-!^ . Moreover Fig. 14 shows that 
the scaled distribution of fh is identical for 3 = and 
5 7^ (the joint probability distributions are however 
very different). It is clear from Figs. 12 - 14, that the 
asymptotic distributions show marked deviations from a 
simple gaussian. To highlight these deviations, we plot 
the scaled log(— log(P(r7ii))) versus log(m^) (Fig. 15), 
a gaussian distribution would have given a straight line 
with slope —1. 




mi/L(t) 

FIG. 12. Scaling plot of the un- normalized 
P(x — mi/L{t)) for 5 = at different times t = 
900(o), 3600(n), 6300(A). Solid line is a fit to Eq. (P). 




mi/L(t) 

FIG. 13. Scaling plot of the un-nornialized P{x — 
mi/L{t)) for q — 0.3 at different times t — 
1350(o), 3600(+), 5400(n). Solid line is a fit to Eq. (|3|). 
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mi / L(t) 

FIG. 14. Scaling plot of the un-normalized P{x — 
mi/L{t,g)) for g = 0(o) and g = 0.3(+) at t = 4500 
showing that the distributions are identical within error 
bars. 



in Figure 12 is an accurate fit to the g ~ data, with 
po = 1, = 1.33 X 10-3 ± 6.0 X 10-'^ p2 = 0.2352 ± 
3.8 X 10-^ P3 = 1.55 X 10-'' ± 1.5 x 10-^ p4 = 5.542 x 
10-3 _^ 7_Q X 10-6. Similarly in Fig. 13, the dark hne is 
an accurate fit to the g — 0.3 data with po = 1; Pi = 
3.95 X 10-3 ± 5.5 X 10-^ = 0.2899 ± 1.3 x 10-^ = 
5.35 xlO-4± 1.3 xlO-^p4 = 1.1913 X 10-2±7.0x 10-6. 
Indeed the odd coefficients are zero to within numerical 
accuracy, indicating that the distribution is even. 

It is conceivable that such corrections would be able 
to salvage the inconsistency issue, since an additive term 
to the right hand side of Eq. ( ^0| ) would not allow us to 
assert that 7/c should be negative for some values of k. 

Though the remedy suggested cures the inconsistency 
problem, it will still give a zero value to the torque con- 
tribution in Eq. (psj), as long as the probability distribu- 
tion of each component of fh is even and independent. 
We have already demonstrated that the single point dis- 
tribution is even, now we shall show that each cartesian 
component of m is independently distributed. 

We numerically calculate P(toi(1), m2(2)) (which we 
label P{x,y)) at equal times ti ^ t2 — t and arbitrary 
separation, say |ri — ra] 4^3 for g = 0.3 (Fig. 16). 
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"-0.2 0,2 0.4 0.6 0,8 1 1.2 1.4 1.6 1.8 



log((m/L(t))2) 

FIG. 15. Deviation of P{x = mi/L) from gaussian 
(straight line) for g = Q. Data have been collected at 
times t = 900(o), 3600(+), 6300(0). 

Figures 12 - 14 suggest that the deviations from gaus- 
sian can be computed by expanding P{m) in a Hermite 
polynomial basis (a strategy advocated in Ref. for 
the scalar (Ising) dynamics), 

oo 

P(x) = ^p„i?„(a;)e-^% (34) 

n=0 

where x = mi(r, i)/\/ S'o(?', t) and Hq{x) — 1, Hi{x) = 
2x and Hn+i{x) — 2xHn{x) — 2nHn-i{x). The dark hne 



'o 




FIG. 16. Normalized joint probability distribution 
P{x,y) where x — toi(1), y = 7712(2) for g — 0.3 at 
t = 2250 and |ri — ra] = 4^3 (averaged over 18 initial 
configurations) . 

To show that the joint distribution is independent in each 
component, we plot the difference /S.{x,y) = P{x,y) — 
P{x)P{y) for g = 0.3 (Fig. 17) and find it to be zero 
within the accuracy of our numerical computation. 
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FIG. 17. Plot of A(a;, y) where x = mi(l), y = 7713(2) 
at t = 2250 and |ri - ra] = 4^3 for g = 0.3. The maxi- 
mum magnitude of A is of the order of errors in A(x, y). 



0.0008 
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0.0004 

y 

0.0002 




-0.0002 



5 10 15 20 25 30 

r 

FIG. 18. y = C^;^{r) (□) and y = C^-ir){*) at 
t = 3600 and r = |ri — for g — are zero within 
the error bars (averaged over 5 initial configurations). 



We are thus forced to admit the second possibility, 
namely that the order parameter (f) cannot be written 
in terms of m alone. For in transforming the spins 4> ex- 
clusively to TO we have implicitly ignored spin waves. A 
most direct demonstration of this is to compare C^^ — 

((^(1) • {(j){2) X V20(2))) with the defect-only contribution 
^3m ~ ("^(1) ■ ("^(2) X V2to(2))) (where rh is computed 
by inverting Eq. (pSh). 

We find that for q = both C„r and C„^ are zero 

^ 30 3m 

within error bars (Fig. 18). This is true even at very early 
times which implies that in the absence of the torque 
the spinwaves decay very fast compared to the relaxation 
timescale of the defects. On the other hand, when g ^ 0, 
we find that the two correlators behave very differently. 
Figure 19 clearly shows that even at late times, C^^ is 
non zero while the defect-only contribution C,^ is zero 

3m 

within errorbars. 

This suggests the following decomposition in terms 
of defect fields (singular part) and spin- waves (smooth 
part), 4> = i7(to) -I- M, when <? ^ 0. Such a decomposition 
gives rise to contributions to Cg^ reflecting the interac- 
tion between defects and spin-waves. 



0.0008 




5 10 15 20 25 30 

r 



FIG. 19. y = C,f{r) (□) and y = C^rJ^){*) at 
t = 3600 and r ~ |ri — r2| for g = 0.3 are dis- 
tinctly different (averaged over 5 initial configurations). 
C^^(r)(+), which has contributions from defects alone, 

3m ^ ^ ^ , , . . 

is zero (within errorbars), whereas C^^{r), which in ad- 
dition involves spin- wave excitations, is non zero. 

We conclude this long meandering section by recount- 
ing its salient results. When g 0, typical spin con- 
figurations at late times consist of slowly moving defects 
and long-lived spin- waves which interact with each other. 
The asymptotic spin distribution cannot be written in 
terms of the distribution of defects alone. When .g = 0, 
the spin waves decay faster, leading to a decoupling of 
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the spin waves and defects at late times. 



with the MSR action written as 



IV. ORDERING DYNAMICS AT T = Tc 

We end this study with a brief discussion of the order- 
ing dynamics Eq. ([^) of Heisenberg spins quenched to the 
critical point. The critical dynamics of this model (called 
Model J in this context) was investigated some time ago 
by Ma and Mazenko On the other hand, the dynami- 
cal renormalization group formalism for quench dynamics 
set up by Janssen et.al., has been used to study Models A 
- C [^Ij] ■ In this section we use the dynamical renormal- 
ization technique to study the quench dynamics of Model 
J (given by Eq. (|l|) at the critical point). Though this 
section does not contain anything new of a fundamental 
nature, it does however compute exponents to all orders 
in perturbation. 

We first demonstrate that the precession term is rele- 
vant for the quench dynamics to Tc- We will then cal- 
culate the z and A exponents at this new fixed point. 
We will show, to all orders in perturbation, that A is ex- 
actly equal to the spatial dimension d. This latter fact, a 
consequence of the conservation law (and indeed true for 
Model B dynamics too), may also be arrived at by the 
general arguments presented in Ref. 

In the absence of the torque term, the nontrivial fixed 
point is given by the Wilson-Fisher value (WF), u* — 
(8/ll)7r^e, where e = 4 — d. Power counting shows that 
the scaling dimension oi g is d/2 + 1 — z + ri/2, where 
the exponents take their WF values z = 4 — r/ and rj = 
(5/242) e^. This implies that the torque g is relevant at 
the WF fixed point when d < 6 [Q. 

We now have to determine this new torque driven fixed 
point and calculate the dynamical exponents z and A. 
Both these exponents can be obtained readily using gen- 
eral arguments, which we briefly discuss. At the new 
fixed point it is clear that g does not get renormalized, 
which implies that z = {d + 2 + ri)/2. Thus a calculation 
of z within perturbation theory reduces to a calculation 
of r] at this fixed point . Likewise A can be obtained 
from the general arguments outlined in Ref. A cru- 
cial ingredient in this argument (valid only for quenches 
to Tc) is the demonstration that S{k, t) obeys a scaling 
form at A: = 0, a feature that was proved in Ref. to 
all orders in perturbation for Model B dynamics. Here 
we directly calculate both z and A using diagrammatic 
perturbation theory, and show that A = d to all orders 
in perturbation. 

This is done within the Martin-Siggia-Rose (MSR) for- 
malism For our problem, the MSR generating func- 
tional is. 



J[(/),0] 



dt / dk<^ 0k ■ 



Ffc 



dki I ^(A:^ - (k - ki)2)4^ X 4_ki 



'k • V^-k 



(36) 



In the expression for the generating functional, the initial 
distribution of the order parameter (gaussian with the 
width = T(7^) enters the form of i/o = / rfk^((/)k(0) • 
'?-k(0)) §. 

Power counting reveals the presence of two different 
upper critical dimensions coming from the quatric term 
(d" = 4) and the cubic torque term (dg = 6) in the action 
J. This implies we have to evaluate the fixed points and 
exponents in a double power series expansion in e = 4 — d 
and e = 6 - d [|. 

The unperturbed correlation C^{ti,t2) and response 
^^•(^11^2) functions, and the bare u and g vertices are 
shown in Fig. 20. Again power counting shows that at 
d = 3, our perturbation expansion does not generate ad- 
ditional terms other than those already contained in J, 
i.e. the theory is renormalizable. However the pertur- 
bation theory gives rise to ultraviolet divergences which 
can be removed by adding counter-terms to the action. 

To remove these divergences, we introduce renormal- 
ization factors (superscripts R and B denote renor- 

malized and bare quantities respectively), 0k (0) = 



TO ^0 



and T^^ = Z. 

Since the dynamics obeys detailed balance, the renor- 
malization factors Z and Z^ are the same as in statics. 
Further the conservation of the order parameter forces 
ZZ = 1 to all orders. 

The new fixed point is given by the zeroes of the (3 
functions of the theory. The (3 functions, calculated from 
the Z factors, get contributions from all diagrams con- 
taining the primitively divergent diagrams F , ~ , F , , - and 

F(,*^, (Fig. 20). 



The new fixed point, to one loop, is given by — 
±Vl927r3e -1-0(63/2)^ u* = {%/l\)TT'^e + 0{€^) (note u* 
does not change from its WF value to all loops) and the 
dynamical exponent z = A — e /2 + 0[e^) 0. 



Z[h,h] = j V{^)V{^) exp I - J[0,0] -i/o[0o] 
+ f dt f dk (4 • l-k + /^k • 0-k) I (35) 
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(a) 



(b) 



C|,(ti,t<i) 



E e(ti-t«)e 



1 , -rk'iti- 12| 2 -rklti+t) 
pie -(T,-i)e ] 




irg£„l^(k, - (k- k,)") - |ru(5a|iS,5 + 5„5 V + 8.75|j8) 



(d) 



(3)_ /Jw _ 



FIG. 20. Unperturbed (a) response function G^, (b) 
correlation function C^, and the (c) two bare vertices 

u and g. Wavy and straight lines represent the 4'u{i) 
and (/'k(i) fields respectively, (d) Primitively divergent 



r(3) andr(^) 



diagrams F^^j , 



The A exponent can be computed from the response 
function Gk(i, 0) = {(f>]^(0) ■ (/)_k(i)) since this is equal to 
the autocorrelation function ro^^(0k(^) • '/'-k(O)), as can 
be seen from the first term in J on integrating by parts. 
The response function gets renormalized by 



G^{t,0) = Z-'/'G^{t,0). 



(37) 



The divergent contributions to Gb could come from 

two sources. Each term in the double perturbation se- 
ries could contain the primitively divergent subdiagrams 
p(2)^ p(3) pj. p(4)^ which we have already accounted for 
by replacing these by their renormalized counterparts. 
The other divergent contribution could arise from the 
primitive divergences of the 1-particlc reducible vertex 
function F(2)(k,t,0), defined by Gk(i,0) = / Gk(i - 
t') F(2) (k, t', 0) dt'. The superficial diver gence of the dia- 
grams contributing to Gk{t, 0) is D = 14(fi — 4) + ■^(d — 
6) — 2 (where Vu (Vg) is the number of u (g) vertices 
respectively). This is negative for all d. For (a) when 
d > 6, the only stable fixed point is the gaussian fixed 
point and so D = —2, (b) when 4 < d < 6, u is irrelevant 
and so D = ^{d - 6) - 2 < and (c) when d < A, D 
is clearly negative. This implies that G^ (t, 0) does not 
get renormalized and Zq ~ I. Consequently A stays at 
its mean-field value of d for this conserved Heisenberg 
dynamics both with and without the torque. 



V. CONCLUSIONS 

Traditional analysis of the asymptotic ordering dynam- 
ics of vector order parameters fociiscs on the dynam- 
ics of defects, and ignores the bulk excitations like spin 
waves which most often decay faster. In this work we 
have looked at the very realistic model of Heisenberg 
spins with precessional dynamics and have shown that 
the longer lived spin waves couple to the defects even at 
late times, driving the system to a new fixed point. The 
new 'torque-driven' fixed point characterized by 2; = 2 
and A « 5.05, is accessed after a crossover time tc ^ 
(where g is the strength of the torque). Crossover scal- 
ing forms describe physical quantities like domain size 
and equal/unequal time correlation functions for all val- 
ues of g. In the absence of the torque, the spin-waves 
decay faster and so do not contribute to the asymptotic 
dynamics. 

We also studied the effects of the torque on the dy- 
namics following a quench to the critical point Tc. The 
torque is relevant with exponents 2; = 4 — e/2 and X = d 
(where e = 6 — d). We found to all orders in perturbation 
theory that \ = d which follows as a consequence of the 
conservation of total magnetization. 

We hope we have provided strong evidence that in or- 
der to go beyond the present approximate theories of 
the asymptotic dynamics of conserved order parameters, 
we need to systematically evaluate contributions coming 
from the interaction of defects with spin-waves. 

We thank D. Dhar for interesting discussions and Y. 
Hatwalne for a critical reading of the manuscript. 



' email : jayajit@rri.ornet.in 
email: madan@rri.crnet.in 

On leave of absence from : Institute of Mathematical Sci- 
ences, CIT Campus, Taramani, Chennai 600 113, India. 

[1] A. J. Bray, Adv. Phys. 43, 357 (1994). 

[2] G. F. Mazenko, Phys. Rev. B 43 , 5747 (1991) ; A. J. 
Bray and K. Humayun, J. Phys. A 25, 2191 (1992) ; 

[3] G. F. Mazenko, Phys. Rev. E 49, 3717 (1994) . 

[4] S. K. Ma and G. F. Mazenko, Phys. Rev. B 11, 4077 
(1975) ; P. C. Hohenberg and B. I. Halperin, Rev. Mod. 
Phys. 49, 436 (1977). 

[5] S. Majumdar, D. Huse and B. Lubachevsky, Phys. Rev. 
Lett. (1995). 

[6] H. K. Janssen, B. Schaub, and B. Schmittmann, Z. Phys. 
B73, 539 (1989) ; H. K. Janssen in From Phase Transi- 
tions to Chaos, eds. G. Gyorgyi et. al. , World Scientific, 
Singapore (1992). 

[7] J. G. Kissnor, Phys. Rev. B 46, 2676 (1992). 

[8] The effects of the torque on the nonconserved dynamics 
of the Heisenberg model was studied in J. Das and M. 
Rao, Phys. Rev. E 57, 5069 (1998). 

[9] M. Seigert and M. Rao, Phys. Rev. Lett. 70, 1956 (1993). 



13 



[10] C. Young, M. Rao and R. Desai, Phys. Rev. E 53, 1 

(1996). 

[11] F. ,J. Wcgner, Phys. Rev. B 5, 4529 (1972) ; C. J. Cama- 
cho and M. E. Fisher, J. Chem. Phys. 94, 5693 (1991). 

[12] D. Stauffer, M. Ferer and M. Wortis, Phys. Rev. Lett. 
29, 345 (1972). 

[13] C. Yeung, A. Shinozaki and Y. Oono, Phys. Rev. E 49, 
2693 (1994). 

[14] M. Abramovitz and I. A. Stegun, Handbook of Mathe- 
matical Functions, (Dover, NY, 1970). 
[15] F. R. Iniguez and A. J. Bray, Phys. Rev. E 51, 188 (1995). 



14 



